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Small-scale structure of the galactic cirrus emission 
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Abstract We examined the Fourier power spectrum characteristics of cirrus structures in 13 sky fields with faint to bright cirrus 
emission observed with ISOPHOT in the 90-200/um wavelength range in order to study variations of the spectral index a. We 
found that a varies from field to field with -5.3 <a <-2.1. It depends on the absolute surface brightness and on the hydrogen 
column density. We also found different spectral indices for the same sky region at different wavelengths. Longer wavelength 
measurements show steeper power spectra. This can be explained by the presence of dust at various temperatures, in particular 
of a cold extended component. For the faintest areas of the far-infrared sky we derive a wavelength independent spectral index 
of a = -2.3+0.6 for the cirrus power spectrum. The application of the correct spectral index is a precondition for the proper 
disentanglement of the cirrus foreground component of the Cosmic Far-Infrared Background and its fluctuations. 
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1. Introduction 

In the far-infrared the most significant component superim- 
posed on the Extragalactic Background is the cirrus emission 
of the Galaxy originating from irregularly shaped interstel- 
lar clouds (Low et al. 1984). The correct determination of 



the spatial frequency fo ( jautier et al., 1992). Since cirrus flue 



the Cosmic Far-InfraRed Background (CFIRB) requires a pre- 



tuations are different from CFIRB fluctuations in at least two 
features (spatial frequency- and brightness-dependence), it is 
possible to separate them from each other in two independent 

ways: 

(1) The decomposition through the brightness dependence 



rise r emoval of this foreground component ( puiderdoni et dj , at a fixed spatial frequency was performed by Kiss et al. (2001 



1997). The CFIRB is characterized by two parameters: the ab- 



solute level and the amplitude of its fluctuations. In recent years 
observations made by ISOPHOT ( f^emke et al. , 1996J), th e pho- 



tometer on board the ISO satellite (Kessler et al., 1996) made 



the determination of the CFIRB fluctuations possible. 

CFIRB fluctua tions are expected to show the following 
properties (see e.g. Hauser & Dwek , 2001 ): (1) the spatial dis- 
tribution corresponds to a flat power spectrum, (2) they are 
isotropic and (3) they have a positive constant value. In con- 
trast, the fluctuations due to galactic cirrus emission scale up 
with absolute brightness and can be described by a multifractal 
structure, yielding a steep Fourier power spectrum. It is usually 
described by one main parameter, the spectral index a: 



hereafter Paper I), and resulted in a new determination of the 
CFIRB fluctuation amplitudes at 90 and 170yimi. 

(2) The decom position using the Fourie r power spectrum 
was introduced by Guiderdoni et al. ( 1997 ). For faint fields, 
where CFIRB fluctuations are expected to be dominant in the 
highest spatial frequency range, the knowledge of the cirrus 



spectral index is mandatory for decomposition. Gautier et al 



p = p x (f /f r 



( 1992 ) determined a spectral index of the galactic cirrus emis- 
sion ff*-3 from the analysis of 100 /mi IR AS scans. For re- 
gions of fa int total FIR brightness emission Lagache & Pugel 
( |200C| ) and |Matsuhara et al.| ( ^000| ) derived CFIRB fluctuation 
amplitudes applying this spectral index of a « -3 for the cirrus 
decomposition. 



(1) C 



where P is the Fourier power, f is the absolute value of the two- 
dimensional spatial frequency and Pq is the Fourier power at 

Send offprint requests to: Cs. Kiss, pkisscs@konkoly.hu 

* Based on observations with ISO, an ESA project with instruments 
funded by ESA Member States (especially the PI countries: France, 
Germany, the Netherlands and the United Kingdom) and with the par- 
ticipation of ISAS and NASA. 



Herbstmeier et al. (1998) analysed a few fields observed by 
ISOPHOT in the 90-180/mi wavelength range. They found a 
variation of spectral indices in the range of -0.5 > a > -3.6, de- 
pending on wavelength, brightness and sky position. However, 
this sample was too small (four fields, two of them measured 
only at one wavelength) for general conclusions, but their re- 
sults clearly indicated that the a — -3 cirrus spectral index may 
not be universal. 
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Figure 1. 100 and 200 yum range 
images of two sample fields in 
t he satellite c o ordin ate system 
( ^aureijs et al. , 2002 ). The cen- 
tral sky positions of the images 
are listed in Table 1. The orien- 
tation of the celestial coordinate 
system is indicated, (a) Cham-II, 
200 fjm (b) Cham-II, lOOyum (c) 
NPC 1, 200 jum (d) NPC 1, 90 /im. 
While the FIR emission in NPC 1 
arises from the same cirrus struc- 
ture at both wavelengths, there 
seems to be an additional com- 
ponent of cold condensations in 
Cham II, as discussed in Sect. 5.1 



Despite the strong impact of the cirrus power spectrum on 
the correct determination of the characteristics of the CFIRB, 
cirrus power spectrum properties were not yet analysed in de- 
tail for A > 100/im and at a resolution achievable by ISOPHOT. 
In this paper we analyse 1 3 sky regions measured by ISOPHOT 
in the 90-200 //m range and investigate the characteristics of 
their Fourier power spectra. 

2. Observations and data analysis 

2.1. Selection of ISOPHOT maps 

We selected 20 FIR maps from the ISO archive, covering 13 
separate sky regions. All selected maps were obtained in the 
PHT22 staring raster observing mode (Laureijs et al., 2002). 



We excluded maps containing strong point sources or obvi- 
ous spatial structures. This selection is a subsample of the 
fields used in Paper I. The selected maps had to show a clear 
sign of cirrus emission in the power spectrum. This requires a 
minimum average surface brightness of (B) » 3 MJysr -1 . Such 
fields are slightly brighter than the faintest fields in the far- 
infrared sky ((B) w 2 MJysr~ ) in the 90-200 //m wavelength 
range after the subtraction of the Zodiacal Light. 

Our sample covers the surface brightness range ~ 3- 
60 MJysr 1 , ranging from faint to bright cirrus, and faint molec- 
ular cloud emission. The basic properties of the maps are sum- 
marized in Table 1. Images of two fields are shown i n Figjj], 
more images are presented in Herbstmeier et al. (1998). 



2.2. Data reduction 

We used the final maps produced for the confusion noise anal- 
ysis in Paper I. For a detailed description of the data reduction 
steps we refer to this paper. The data reduction comprised the 
following main steps: 



- Basic data analysis with PIA[] V8.2 ( pabriel et ajj , |1997[ ) 
from ERD to AAP level 

- Flat-fielding using first quartile normalization 

- Subtraction of the Zodiacal emission 

2.3. Fourier power spectrum 

Fourier transformation of the maps was performed by rou- 
tines written in IDL^| based on th e standard IDL FFT routine 



with periodogram normalisation (Press et al., 1992). For all 



maps the highest spatial sampling frequency corresponds to the 
pixel size of the CI 00 and C200 array cameras. The Nyquist 
limit which allows correct sampling (without aliasing) is half 
of that frequency ( foess et~aT , 1992, Chapter 12.1), i.e. the 
double pixel size of the array. Therefore, the upper limit of 
the spatial frequency considered for the Fourier power spec- 



PIA is a joint develpoment by the ESA Astrophysics Division 
and the ISOPHOT consortium led by the Max-Planck-Institut fur 
Astronomie (MPIA), Heidelberg 

2 Interactive Data Language, Version 5.2 and 5.3, Research Systems 
Inc. 
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NGP 


180 


46.0x46.0 


88.9 


73.0 


3.1+0.1 


-2.53±0.22 


1.8x10 s 


1.59 


0.95 
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9.0xl0 19 


North Galactic Pole 


Juvela et al. (2000) 




























Herbstmeier et al. (1998) 


Draco 


90 


7.7x7.7 


89.8 


38.6 


5.7±0.3 


-2.45+0.17 


5.3xl0 3 


1.47 


0.76 


- 


2.1xl0 20 


Draco nebula 


Herbstmeier et al. (1998) 


Draco 


170 


27.6x18.4 


89.8 


38.6 


4.7±0.5 


-4.09+0.22 


2.9xl0 3 


2.71 


1.39 


- 


2.1X10 20 


Draco nebula 


Herbstmeier et al. (1998) 


Cep 


90 


14.7x20.7 


108.0 


13.6 


23.5+3.1 


-3.42+0.15 


3.3xl0 4 


1.96 


1.01 


6.9 


2.5xl0 21 


Cepheus Flare 


Herbstmeier et al. (1998) 


Cep 


170 


21.5x24.5 


108.0 


13.6 


68.2±12.7 


-4.19+0.12 


3.0xl0 5 


1.95 


1.56 


6.9 


2.5xl0 21 


Cepheus Flare 


Herbstmeier et al. (1998) 


NPC1 


90 


29.9x29.1 


121.6 


24.2 


7.6±0.8 


-3.10+0.04 


5.0xl0 3 


2.16 


0.64 


6.6 


6.9xl0 20 


Polaris Flare 


Gautieretal. (1992) 


NPC1 


200 


30.7x29.1 


121.6 


24.2 


26.1+5.3 


-3.46+0.09 


1.1x10 s 


2.30 


0.95 


6.6 


6.9x1 20 


Polaris Flare 


Gautier et al. (1992) 


NPC2 


90 


29.9x29.1 


122.0 


24.6 


7.4±0.7 


-2.83+0.09 


5.7xl0 3 


1.74 


0.89 


5.5 


6.7xl0 20 


Polaris Flare 


Gautier et al. (1992) 


NPC2 


200 


30.7x29.1 


122.0 


24.6 


22.6±3.8 


-3.39+0.08 


9.7xl0 4 


2.19 


0.63 


5.5 


6.7xl0 20 


Polaris Flare 


Gautier et al. (1992) 


North 1 


90 


29.9x29.9 


100.0 


14.8 


13.3+1.4 


-3.05+0.09 


9.7xl0 3 


1.83 


0.91 


4.5 


1.6xl0 21 


LDN1122 


Yonekura et al. (1997) 


North 1 


200 


30.7x30.7 


100.0 


14.8 


34.4±5.6 


^1.70+0.17 


4.5xl0 4 


2.43 


0.52 


4.5 


1.6xl0 21 


LDN1122 


Yonekura et al. (1997) 


North2 


90 


29.9x29.9 


108.0 


15.2 


12.4+1.8 


-3.65+0.13 


9.2xl0 3 


1.42 


0.90 


6.6 


2.1xl0 21 


LDN1147, 1148 


Lee et al. (2001) 




























Yonekura et al. (1997) 


North2 


200 


30.7x30.7 


108.0 


15.2 


40.7±12.5 


-3.81+0.09 


3.2xl0 5 


1.94 


0.61 


6.6 


2.1xl0 21 


LDN1147, 1148 


Lee et al. (2001) 




























Yonekura et al. (1997) 


M01 


180 


27.6x27.6 


100.0 


30.6 


5.0±0.2 


-2.56+0.41 


2.3xl0 2 


0.89 


0.75 


- 


4.2xl0 20 


faint field in Dra 


Herbstmeier et al. (1998) 




























Abraham et al. (1997) 


M03 


180 


27.6x27.6 


117.6 


46.1 


3.7±0.1 


-2.12+0.23 


5.8xl0 2 


0.82 


0.71 




1.5xl0 20 


faint field in UMi 


Abraham et al. (1997) 


TMC2-1 


200 


30.7x30.7 


173.9 


-15.7 


56.5±8.9 


-5.26+0.15 


3.7xl0 5 


1.74 


0.86 


18.9 


1.7xl0 21 


LDN1529, 1531 


Lee et al. (2001) 


























Taurus 2 mol. cl. 


Chappell & Scalo (2001) 


TMC2-2 


200 


30.7x30.7 


174.3 


-15.9 


63.7±8.5 


-3.49+0.14 


3.8xl0 5 


1.66 


1.04 


16.3 


1.8xl0 21 


LDN1529, 1531 


Lee etal. (2001) 


























Taurus 2 mol. cl. 


Chappell & Scalo (2001) 


ChalS 


100 


19.9x19.9 


297.3 


-16.2 


29.3±3.2 


-3.80+0.08 


1.4xl0 5 


2.47 


0.79 


10.3 


9.2xl0 20 


Chamaeleon mol. cl. 


Toth et al. (2000) 


Chamll 


100 


31.4x31.4 


303.5 


-14.2 


23.9±3.3 


-2.95+0.08 


1.7x10 s 


0.91 


0.89 


16.1 


1.2xl0 21 


DCld 303.3-14.3, 


Toth et al. (2000) 


Chamll 


200 


30.7x30.7 


303.5 


-14.2 


38.4±14.6 


-5.17+0.18 


2.4xl0 5 


2.48 


0.53 


16.1 


1.2xl0 21 


303.5-14.4 


Toth et al. (2000) 



Table 1. Basic observational and derived parameters of fields mapped by ISOPHOT. The columns are: (1) name of the field as specified by the observer; (2) central wavelength of filter; 
(3) map size; (4)-(5) center of the field in galactic coordinates; (6) median surface brightness from ISOPHOT photometry with Zodiacal Light subtracted; quoted standard deviations 
are a measure of the dynamic range of the map. (7) spectral index of the power spectrum; (8) fluctuation power at the reference spatial frequency 1/do (do = 4'); (9)-(10) parameters 
describing the relative strength of the instrument noise power spectrum to the signal power spectrum (see Eq. 2 for the definitions); (11) average neutral hydrogen column density 
(Dickey & Lockman, 1990); (12) 12 CO integrated intensity (Dame et al., 2001); (13) associated objects or dark clouds located (at least partially) inside the field. (14) references with 
the description of the field. 
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trum was set to the Nyquist-limit, which is min = 92" for im- 
ages with A < lOOyum (CI 00 detector) and 6 min = 184" for im- 
ages with A> 100 yum (C200 detector). Since the pixel sizes of 
the ISOPHOT CI 00 and C200 cameras were adjusted to the 
diffraction limit, this limit should be kept, even if the raster 
sampling were finer than the pixel size. 

Instead of averaging P(f) fluctuation power values in an- 
nuli for each f = |f| spatial frequency, we used all individual 
data points (f, P(f) pairs) to derive the final power spectrum. 
The spectral index a is derived by robust line fitting to all data 
points in the log(f)-log(P) space. 

In most cases the whole power spectrum could be satisfac- 
torily fitted by one spectral index. For some faint fields, how- 
ever, the high frequency end of the spectrum showed a different 
(flatter) spectral index, due to a higher structure noise level. In 
such cases we only used the low and mid frequency part to de- 
rive a. The reason for this flattening will be discussed below in 
Sect. 4. 
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3. Features affecting the power spectrum 

There are several contributions to the power spectrum that can 
affect its final shape. We analyse the effects of the instrument 
noise, a point-source footprint and an extended source pro- 
file on the final power spectrum. The correct separation of 
the power spectrum components requires well-distinguishable 
spectral indices, otherwise they cannot be disentangled. 

3.1. Power spectrum of the instrument noise 

Instrument noise can be a significant contribution to the fluctu- 
ation power at low surface brightness and high spatial frequen- 
cies. As discussed in detail in Paper I, the most representative 
estimates of the instrument noise are PIA-noise (which reflects 
the statistical uncertainties in signal derivation) and flat-field 
noise (which reflects the variation of the signal from individual 
detector pixels at the same sky position). Flat-field noise can 
only be determined for oversampled maps. 

We assume that both the signal and instrument noise power 
spectra - derived from the surface brightness maps and the cor- 
responding noise maps, respectively - are described by a power 
law behaviour as in Eq. [l]. Here the fo values correspond to the 
Nyquist-limits, i.e6» ml „ = 92" and 184" for the C100 and C200 
detectors, respectively. 

The signal power spectrum is not affected by the instru- 
ment noise, if the power spectrum of the instrument noise is 
reasonably weaker than that of the signal in the frequency range 
analysed. This is ensured, if two criteria are fulfilled simulta- 
neously: (1) the fluctuation power of the instrument noise P, is 
lower by at least a factor of 2 than that of the signal Pi at the 
resolution limit, 6> m ; n and (2) the power spectrum of the signal is 
steeper than that of the instrument noise. These conditions are 
parametrized by the following two quantities (i and s denote 
instrument noise and signal, respectively): 



0_ 1 



-1 




6a — a\ - a s > £Po 



logio(^) 



> 0.3 



(2) 



1 

6a 



Figure 2. Relative strength of the instrument noise- to the sig- 
nal power spectrum, (a) Subsample of 56 maps from Paper I., 
observed with the C200 detector in oversampling mode. Black 
dots: flat-field noise; gray dots: PIA-noise. (b) Subsample of 
60 maps from Paper I., observed with the C100 detector. In this 
case no representative values of the flat-field noise could be 
dervied due to the lack of a sufficient number of oversampled 
maps. The average values and the standard deviations both in 
5a and 5Po are marked via large dots and error bars, respec- 
tively. 



For the instrument noise analysis we used a larger sample 
of maps than for the cirrus power spectrum study. In Table 1 
we provide the respective values of 6a and 6Pq for our special 
fields, which can be compared with the ranges of 5a and 5Po 
for the larger sample of maps. 

In the case of the C200 detector we calculated 6a and 6Po 
for a subsample of 56 maps from Paper I, observed with full 
oversampling (Fig. ||a). 
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Most of the flat-field and PIA-noise data points are located 
in the quadrant where both 6a and 6Po are positive. Although 
there are a few points in the negative 6a regime, they are close 
to zero, representing nearly parallel instrument noise and signal 
power spectra. Since in all cases 6Pq > 0.3, the instrument noise 
power spectrum does not exceed the signal power spectrum for 
frequencies below 6~\ . The data points with a S= belong to 
the faintest regions of the sky with a nearly flat power spectrum, 
showing no signs of cirrus (a s and c,- * 0). The distributions 
of the PIA- and flat-field noise data points are rather similar, 
and the average values are within the standard deviations in 5a 
and 6Pq (big black and gray dots with error bars in Fig. ||a). 
This demonstrates that the method of noise determination has 
no impact on the final result. 

In the case of the C100 detector it was not possible to per- 
form a similar investigation due to the lack of suitable over- 
sampled maps. Therefore, we analysed 60 observations from 
Paper I by calculating the PIA-noise values only (Fig. ||b). Also 
the CI 00 detector measurements fulfil the requirements of a 
suitable signal-to-noise ratio. 

This noise analysis justifies, that for both C200 and CI 00 
maps of Table 1 instrument noise does not affect the final 
shape of the power spectrum, in particular since for all of them 
6a>0.8 and 5P > 0.5. 



3.2. Power spectra of point- and extended sources 

In Fig. [|a we present a model cirrus spectrum with an expected 
cirrus spectral index of a = -3 convolved with the 90 fim 
and 200 yum theoretical footprints (footprint = convolution of 
the telescope point spread function with the pixel aperture, 



Laureijs et al. , 2002 ). This figure shows, that the footprint has 
only a little effect below the actual Nyquist-limit, and this ef- 
fect can be neglected when deriving the spectral index. A strong 
point source may dominate the whole spectrum, but maps con- 
taining these features were not included in our sample. 

As will be shown in Sect. 5.1, extended sources with ex- 
cess emission may be superimposed on the cirrus-like intensity 
distribution. We tested the effect of these sources on the power 
spectrum assuming a Gaussian shape, and a FWHM of 1/3 and 
1/6 of the size of the original map, respectively, as presented in 
Fig. ||b. The shape of the resulting power spectrum strongly 
depends on the relative strength of the cirrus and extended 
source components. The typical effect of the extended source 
is the appearance of a 'hump' at medium spatial frequencies, 
if the strengths of the two components are comparable. Similar 
(however weaker) features can be identified in some of our real 
power spectra, especially at 200 /im (see Fig. The high fre- 
quency part of an extended source power spectrum is rather 
steep and therefore - if dominant - can lead to a power spec- 
trum steeper than the one of pure cirrus. A set of extended 
sources has a power spectrum similar to that of a single one, 
but the higher the number of sources, the shallower the final 
power spectrum. In certain configurations these spectra may be 
described by an average spectral index of a « -3, leading to a 
spectrum, which is in practice undistinguishable from the one 
of cirrus. 
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Figure 3. Schematic representation of the contribution of dif- 
ferent emission components to the final power spectrum, 
(a) Convolution of a theoretical cirrus spectrum of ar- 
bitrary stength (a = -3, P=10 3 Jy 2 sr -1 at f = 0.1 arcmin -1 ) 
with ISOPHOT's theoretical point-source footprint at 90 fim 
(dashed line) and at 200 yt/m (dash-dotted line). The pure, 
non-convolved cirrus spectrum is also shown (solid line), (b) 
Resulting power spectrum of a cirrus-like intensity distribution 
(a = -3, solid line) and an extended source of Gaussian shape 
with a FWHM of 1/6 (dotted line) and 1/3 (dash-triple-dotted 
line) of the size of the original map. The Gaussians are normal- 
ized to P = 10 5 Jy 2 sr _I atf=0 arcmin . The Nyquist-limits for 
the CI 00 and C200 filter maps are marked both in (a) and (b) 
by vertical dashed and dash-dotted lines, respectively. 



4. Results 

Table 1 and Fig. ^ present the results of the power spectrum 
analysis. 

The fields in our analysis were characterized by the follow- 
ing physical parameters: median FIR surface brightness from 
ISOPHOT photometry with Zodiacal L ight subtracted (this 
work ), neutral hydrogen column dens ity (Dickey & Lockman , 
19901 ) and 12 CO integrated intensity (|Dame et al.|, |200l|). The 
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Figure 4. Power spectra of six sample fields, two of them observed by both the CI 00 and the C200 camera. The images cor- 
responding to the first four power spectra are presented in Fig. [j]. Images of M01 and NGP can be found in Herbstmeier et al 



( ]1998| ). (a) Cham-n, 200 jum (b) Cham-II, 100 jum (c) NPC 1, 200 jum (d) NPC 1, 90 //m (e) M01, 180/im (f) NGP, 180/rni. Each 
dot represents a (f,P(f)) pair, the solid line is smoothed average to outline the trend. The high frequency end is determined by the 
Nyquist-limit, the low frequency end by the size of the map. 



power spectra follow quite well the expected power-law be- 
haviour up to the spatial frequency corresponding to the res- 
olution limit, and therefore can be satisfactorily fitted by one 
single spectral index a. The only exceptions are the fields NGP 
(Fig. @f), M01 (Fig. §0 and M03 at 180/rni. They show a flat- 
tening of the spectrum for the highest spatial frequencies which 
was eventually excluded in the derivation of the a values given 
in Table 1 . 

We suggest that the relatively strong structure noise at high 
spatial frequencies arises from the fluctuations due to the extra- 
galactic background. As discussed in Paper I, the extragalac- 
tic background is the major contributor to the structure noise 
at high spatial frequencies for the faintest regions of the far- 
infrared sky measured with the C200 detector. These three 
faint fields have an average surface brightness comparable to 
that of the cosmological fields. The extragalactic background 
is expected to have a fluctuation power of ~ 8xl0 2 Jy 2 sr~ l at 
170 /mi with our FFT normalization. Fluctuation powers at the 
reference spatial frequency listed in Table 1 (column #8) are fit- 
ted values using the cirrus-like power spectrum, therefore they 
can be lower than the expected CFIRB fluctuation level. 

We calculated the median values of the spectral indices for 
the CI 00 and C200 fields («ci and aci, respectively), in or- 



der to compare them with previous results. This resulted in 
act = -3.15+0.48 and a C 2 = -3.87+1.06. The median C100 
spe ctral index is not mu ch different from the one determined 
by |Gautier et"al| ( |1992| ) from the analysis of IRAS 100/rni 
sca ns, ho wever it shows a relatively large scatter. Abergel et 
al. ( 1996| ) performed a pixelwise determination of the spectral 
index for a 12?5xl2?5 ISSA map close to the South Ecliptic 
Pole. At lOOjt/m they obtained an average spectral index of 
c?ioo = -3.34 (using the relation a = -2(J3 - 1) between their 
/3-value and our spectral index). The distribution of the spectral 
index over the field had a Gaussian shape with cr a = 0.36. This 
means, that within one field the variation is comparable to the 
variation we found for our sample fields. The C200 spectral in- 
dices are significantly steeper with an even higher scatter than 
the CI 00 spectral indices. 

Spectral indices show a clear dependence both on the sur- 
face brightness and on the neutral hydrogen column density. 
Brighter fields and fields with higher HI column density usually 
have a steeper power spectrum (Fig. ^]). We tested the degree of 
correlation by calculating the linear Pearson correlation coeffi- 
cient which resulted in 0.69, 0.80 and 0.67 for the logio(Bci)- 
ffci, logio(Bc2)-Q'c2 and log ioN(HI)-a relations, respectively, 
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without the data points of the Draco region (marked in Fig. |] 
by an arrow). 

If we exclude in Fig. 5b the three data points with high 
Bc2 and a < -4.5, then the logio(Bc2)-ac2 correlation is much 
tighter with a linear Pearson correlation coefficient of 0.94. 
This exclusion is justified by the fact, that all three fields have a 
high HI column density and two of them show a large acz/aci 
ratio (Fig. 5c), indicative of a relatively large molecular hydro- 
gen content, as discussed below in Sect. 5.2. 

It is possible to determine an aci value for the faintest fields 
of the far-infrared sky, by a robust linear extrapolation of the 
logio(Bc2)-Q'c2 relation. In order to account for the contribu- 
tion by the CFIRB the Bc2 values were reduced beforehand 
by 0.8MJysr~' (Paper I). For all points except the Draco re- 
gion this resulted in aci = -2.3+0.6 for the faintest cirrus fields 
((B) « 2 MJysr -1 ). Excluding in addition the three points with 
a<-4.5, aci =-2.1+0.4. The impact of this power law in 
separating the CFIRB from the galactic cirrus is discussed in 
Sect. U 

In Fig. ^| we plot the ratio of the spectral indices aci/aci 
determined both for a C200 and a CI 00 map of the same field. 
Assuming that the FIR radiation at 100 and 200 yum is emitted 
by the same structures the ratio aci/aci should be as 1. It is 
obvious, that acilaci varies in the range 1.0-1.8. Since the de- 
termination of the spectral index is quite accurate (see Table 1), 
this deviation cannot be explained by measurement uncertain- 



ties. We propose an explanation in Sect. 5.1 



5. Discussion 

5.1. Wavelength dependence of the spectral index 

Galactic cirrus emission is expected to have a spectral energy 
distribution of a modified black body, with v 2 emissivity law. 
and a dust colour temperature of 17.5±1.5K ( |^agache et al 



1998). If the FIR emission of a certain cloud can be well de- 
scribed by a single dust temperature, then the emission mea- 
sured at two separate wavelengths should be strongly corre- 
lated, showing the same spatial structure and the same power 
spectrum. On the other hand, the presence of another compo- 
nent with different colour temperature can lead to an excess 
intensity at one wavelength. This yields also a different spatial 
structure of the emission, and therefore different power spectra 
at two wavelengths. Indications for such a dichotomy in cloud 



hierarchy were found by Abergel et al. ( 1996 ) from their analy 



sis of the 60 and lOOyimi IRAS maps of the South Ecliptic Pole 
field. Uncorrelated features come from the coldest regions of 
the analysed fields. 

We tested the uniqueness of the colour temperatures for 
fields, where measurements at two wavelengths were avail- 
able. The colour temperature is represented by the flux ratio 
of the two wavelengths. Two examples are presented in Fig. [7|. 
The 90-100 jum surface brightness maps were smoothed to the 
resolution of the 1 70-200 yum measurements in order to pro- 
duce Bci - Bc2 pixel-to-pixel scatter plots, where Bci and Bc2 
are the surface brightnesses of the 90-100 /vm and 1 70-200 /im 
map pixels, respectively. In most cases the scatter plot showed 
a well correlated Bci - Bc2 emission and could be fitted by 
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Figure 5. Relationship of the spectral index with (a) the 100 /-im 
range (90-100yum) median surface brightness Bci, (b) the 
200 /im range (1 70-200 yum) median surface brightness Bc2, 
(c) the average neutral hydrogen column density. Crosses 
and open circles mark fields observed at A < 100 yum and 
170 < A < 200 /im, respectively. The data points of the Draco 
170 /im field are marked by arrows. Error bars in (a) and (b) 
represent the standard deviations in Bci and Bc2 and they mark 
the errors in fitting the spectral index in a. 
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Figure 6. Ratio of the spectral indices of sky regions mea- 
sured in both a long wavelength (1 70-200 fim, aci) and at short 
wavelength (90- 100 /mi, aci) filter versus the average surface 
brightness in the C100 filter (90 or 100/mi, Bci). The regions 
are labelled by the names as given by the observer and as listed 
in Table 1 . 

one straight line, corresponding to one single dust temperature. 
Using this relation we derived an excess surface brightness AB 
map, by subtracting the scaled up Bci emission from Bc2- 



AB = B a -[Ax xB cl +A ] 



(3) 



where Ai and Ao are the slope and the interception of the linear 
fit to the scatter plot, respectively. 

However, as presented in Fig. [7]a, in some cases there was 
more than one characteristic colour temperature present in a 
field. The most conspicuous example is the Cham II field (gray 
dots in Fig. 0a). In this case, after removing the cirrus con- 
tribution according to the procedure of Eq. || the excess sur- 
face brightness image AB still showed the steep index of the 
original map. The steep AB power spectrum can be associ- 
ated to a br ight extended source, as suggested by the examples 
in Sect 3.2 and, indeed, the region contains dark cloud cores, 
namely DCld 303.3-14 and DCld 303.5-14.4 (see Table 1). The 
expected flattening at low spatial frequencies due to the ex- 
tended source power spectrum (see Sect. 3.2) can also be iden- 
tified. The characteristic size of the extended source derived 
from the power spectrum agrees well with the typical cold core 
extent of ~10' found by Toth et al. (2000) in Chamaeleon, with 
temperatures of < 15 K. A similar effect can be observed for 
the North 1 field, however, the differences are smaller than in 
the case of the Cham II region. The feature in the power spec- 
trum can be identified with the molecular ( 13 CO) core of the 
dark cloud LDN 1122 (lYonekura et all |1997h. In the case of 



the NPC 1 field we found that the power spectrum of the ex- 
cess surface brightness image is significantly shallower than 
that of the 200 /mi surface brightness image, because this lacks 
any dark cloud core. 

No colour temperature could be derived for the Draco re- 
gion via the scatter plot procedure, due to its small spatial ex- 
tent and narrow surface brightness range. However, the pres- 
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Figure 7. Examples to demonstrate the effect of a non-uniform 
color temperature within the same field, (a) C100 vs. C200 sur- 
face brightness scatter plots. The slope is closely related to the 
color temperature (see text for details). Black and gray dots rep- 
resent the data points of the North 1 and Cham II fields, respec- 
tively, (b) Power spectra of the AB images (cirrus component 
subtracted) of Cham II (gray) and North 1 (black). 



ence of a relatively bright extended source i s obvious in the 
170/mi image (see j^erbstmeier et al , 1998 ). This produces 
very likely the very steep power spectrum observed at this 
wavelength, despite the low av erage surface brightness (see 
Fig. ||), as already mentioned by Herbstmeier et al. ( 1998 ). 

These findings suggest that the wavelength dependence 
of the spectral indices is caused by the combined effects of 
non-uniform dust temperatures and the presence of extended 
sources in the fields we selected for our study. The separation of 
the cirrus and extended source components in the power spec- 
tra can be more easily performed by including spatial scales 
larger than the characteristic extent of our maps. At 170/mi 
such an analysis would be possible by utilizing the ISOPHOT 
Serendipity Survey slews with a typical length of a few degrees 
(jBogun et al], |1996|). 
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5.2. The effect of atomic-to-molecular phase transition 

The correlation of the spectral indices with surface bright- 
ness and with neutral hydrogen column density is obvious 
for faint fields, while a scatter in a is observed for bright 
and high HI column density fields, as presented in Fig. H, 
These fields also show strong molecular ( 12 CO) emission, as 
presented in Table 1. This large scatter can be observed for 
N(HI) > 10 21 crrT 2 , which is approximately the column density 
of the atomic-to-molecular phase transition in the galactic cir- 



rus (deVries & van Dishoek, 1988; Meyerdierks et al. 1990) 



In these clouds the ratio of molecular hydrogen column den- 
sity to atomic hydrogen col umn density ma y vary in the range 
of N(H 2 )/N(HI)« 0.2-5.0 (|Magnani et ak|, |1985|). Therefore 



clouds showing the same N(HI) may have different total hy- 
drogen masses. This diverse molecular content - and therefore 
diverse final structure - can lead to the observed scatter in the 
spectral indices. 

5.3. Implication of cirrus power spectra on the CFIRB 
fluctuations 

The spectral index for faint fields observed in the 1 70-200 /im 
wavelength range is greater than t he usually ass u med a = -3, 
which was originally derived by 3autier et al. (1992) from 



100 /im IRAS scans of medium-brightness cirrus fields. Due 
to the lack of a better determination this was widely applied to 
the faintest fields of the far-infrared sky and for longer wave- 
lengths (e.g. 3uiderdoni et al. , 1997; pTgache & Puget , 2000). 
We showed that brighter extended sources - which are often 
present in cirrus fields - may significantly increase the steep- 
ness of the power spectrum. Therefore it is plausible, that the 
a — -3 power law found for medium surface brightness fields 
results from the superposition of the power spectrum of a less 
steeper fractal structure with extended sources. Our analysis 
provides a = -2.3+0.6 for the spectral index of the faintest cir- 
rus fields. Although there is a discrepancy between a-values 
determined at different wavelengths for brighter fields, the 
spectral indices should be the same for faint fields, since these 
lack additional extended sources and multiple dust colour tem- 
peratures. We finally conclude that a uniform a = -2.3+0.6 cir- 
rus spectral index - independent of wavelength - has to be ap- 
plied for the faintest areas ((B) » 2 MJysr -1 ) of the far-infrared 
sky. The application of this a-value would decrease the CFIRB 
fluctuation amplitudes obtained by previous works in the 170— 
200 /mi range by 5-20% (Kiss et al., 2002, in prep.). For some- 
what brighter regions it is mandatory to determine the proper 
spectral index from the properties of the local structures. 

6. Conclusions 



We examined the Fourier power spectrum characteristics of 13 
sky fields with faint to bright cirrus emission in the 90-200/im 
wavelength range and derived the spectral index of the power 
spectrum, a. We found that a varies from field to field. It has a 
clear dependence on the absolute surface brightness and on the 
hydrogen column density of the field. For longer wavelengths 
(A > lOOyum) the scatter in a is larger, especially for higher hy- 



drogen column densities. This is because for N(HI) > 10 21 cm -2 
the atomic-to-molecular phase transition leads to various a 
values especially in the 170-200 yum range, where the diverse 
molecular content is easier to observe. We also found a sig- 
nificant difference in spectral indices for the same sky region at 
different wavelengths. At longer wavelengths (A > 100 fim) the 
power spectra are steeper. This can be explained by the coexis- 
tence of dust components with various temperatures within the 
same field and cold extended emission features which signifi- 
cantly affect the power spectrum. A spectral index of the cirrus 
power spectrum with a value of a = -2.3+0.6 - independent of 
wavelength - could be derived for the faintest areas of the far- 
infrared sky. The precise determination of the spectral index 
for each field will enable to correctly disentangle the galactic 
foreground components of the CFIRB. 
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